% 利用FFT实现快速线性卷积
clear;
% 卷积信号
N1 = 16;
N2 = 17;
n1 = 0:1:N1-1;
n2 = 0:1:N2-1;
xn = (1.0/3).^n1;
hn = 0.5.^n2;

% 自定义FFT滤波器
yn = myfilter(xn, hn);
yn2 = conv(xn,hn);
% 卷积结果长度
n = 0:N1+N2-2;
subplot(2,1,1);
stem(n,yn);
title('FFT快速卷积结果');
subplot(2,1,2);
stem(n,yn2);
title('conv函数卷积结果');


